A study of qq states in transverse lattice QCD using alternative fermion formulations 
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In this work we investigate qq spectra and wavefunctions of light front transverse lattice Hamiltonians that 
result from different methods of formulating fermions on the transverse lattice. We adopt the one link approx- 
imation for the transverse lattice and Discrete Light Cone Quantization to handle longitudinal dynamics. We 
perform a detailed study of the continuum limit of DLCQ and associated techniques to manage severe light front 
infrared divergences. We explore the effects of various parameters of the theory, especially, the strength of the 
£f>^ • helicity-flip interaction and the link mass on spectra and wavefunctions. 
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• I. INTRODUCTION 

<n : 

A promising method to calculate observables in QCD is the transverse lattice formulation OS Ht- In this method, one 
keeps x ± = x° ±x 3 continuous and discretizes the transverse space spanned by coordinates x x = (x l ,x 2 ). With the gauge 
£^ ' choice A + = 0, A~ becomes a constrained variable which can be eliminated in favor of dynamical gauge variables. So far very 

encouraging results have been obtained in the pure gauge and meson sectors 
£T) • Due to the doubling phenomena, fermions on the lattice pose challenging problems. To date, calculations of meson properties 
using transverse lattice have employed Wilson fermions |7[ . It is well-known that the Wilson term explicitly breaks the chiral 
symmetry and makes it difficult to explore the consequences of spontaneous chiral symmetry breaking in the chiral limit. In this 
limit, in one link approximation on the transverse lattice, the Wilson term can be adjusted to produce the desirable level splitting 
between % and p. This, however, results in the undesirable consequence that the splitting of the p multiplet is almost as large 
rSlL as the % — p splitting. Because of the doubling problem, one cannot keep the Wilson term very small. Thus it is desirable to 
explore other formulations of fermions on the transverse lattice that may have different chiral properties. 

In a recent work IS] we have addressed the problems of fermions on the light front transverse lattice. We proposed and 
numerically investigated different approaches of formulating fermions on the transverse lattice. In one approach, which uses 
■ forward and backward derivatives, fermion doubling is absent and the helicity flip term which is proportional to fermion mass in 
light front QCD becomes an irrelevant term in the free field limit. In the literature, symmetric derivatives have been used which 
leads to fermion doubling due to the decoupling of even and odd lattices. Using the light front staggered fermion formulation 
and the Wilson fermion formulation, we studied the removal of doublers from the spectrum. Our investigations lead to the 
identification of an even-odd helicity flip symmetry of the light front transverse Hamiltonian, the absence of which means the 
removal of doublers in all the cases that we studied. 

In this work we make a detailed comparison of various light front QCD Hamiltonians that result from different ways of 
formulating fermions on the transverse lattice. As the first step in our calculations, we adopt the one link approximation in the 
meson sector which has been widely used in the literature. (Only very recently, the effect of additional links in the meson sector 
has been investigated |6]). Since the one link approximation is very crude, our aim is not to reproduce physical observables. 
Rather, we explore the effects of various coupling strengths on the low-lying spectra and wave functions and compare two 
different formulations. 

We use Discretized Light Cone Quantization (DLCQ) |9] to address longitudinal dynamics. Because of the presence of severe 
light front infrared divergences, a major concern here is the reliability of DLCQ results when calculations are done at finite 
resolution K and results are extrapolated to the continuum (K — > °°). In meson calculations so far, K < 20 have been chosen. In 
this work we perform a detailed study of the continuum limit of DLCQ by performing calculations at larger values of K. 
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In the meson sector, in the zero link approximation, at each transverse location we have a two-dimensional field theory which 
in the large N c limit (where N c is the number of colors) is nothing but the 't Hooft model. In this well-studied model, excited 
states are simply excitations of the qq pair, which contain nodes in the wavefunctions. The picture changes when one link is 
included thereby allowing fermions to hop. The admixture of qq link states with qq states is controlled by the strengths of the 
particle number changing interactions and the mass of the link field. One link approximation is a priori justified for very massive 
links and/or weak particle changing interaction since in this case low lying excited states are also qq excitations. Likewise, for 
large particle changing interaction strength and/or light link mass, low lying excited states are qq link states. We explore the 
spectra and wavefunctions resulting from the choice of various regions of parameter space. 

The plan of this paper is as follows. In Sec. II we present the details of the light front transverse lattice Hamiltonian resulting 
from the use of forward and backward derivatives and the resulting effective Hamiltonian when the unitary link variables are 
replaced by general complex matrices. In this section we also present the canonical transverse lattice QCD Hamiltonian resulting 
from the addition of the Wilson term. Sec. Ill contains comparisons of numerical results for the two Hamiltonians. Finally Sec. 
IV contains our summary and conclusions. Typical terms in the Hamiltonian with forward and backward derivatives in the 
Fock representation in DLCQ is presented in Appendix lAl Explicit expressions for the states are given in Appendix 151 For 
completeness, explicit expressions for the matrix elements in the forward-backward case and the Wilson case are presented in 
AppendiceslClandlDl 



II. HAMILTONIANS 



Due to the constraint equation in light front theory, different methods are possible to put fermions on the transverse lattice. In 
this section we present the detailed structure of the resulting QCD Hamiltonians for two methods studied in Ref. |8], namely, 
forward and backward derivatives and symmetric derivatives together with Wilson term. 



A. Hamiltonian with forward and backward derivatives 

Details of the derivation of the fermionic part of the Hamiltonian are already given in Ref. 0. Here we give the details of the 
gauge field part of the Hamiltonian. Non-linear constraints on the unitary link variables make it difficult to perform canonical 
quantization. We also present the effective Hamiltonian when non-linear unitary variables are replaced by linear variables. 



1. Gauge field part of the Lagrangian density 



The gauge field part of the Lagrangian density in the continuum is 



2g z 



(2.1) 



where F pa = d p A a — d a A p + [A p ,A a ] withA? = igA pa T a . Here p, a = 0, 1,2,3 and a = 1,2, ... ,8. For ease of notation we 
suppress the dependence of field variables on the longitudinal coordinate in this section. With the gauge choice A + = 0, the 
Lagrangian density can be separated into three parts, 



Here S£t depends entirely on the lattice gauge field U r (x). 

JCt = -TT Y[ I Tr\u r (x)U s (x + ar)U- r (x + ar + as)U- s (x + as) - 1 
The purely longitudinal part Jz?l depends on the constrained gauge field A - , 



^L = ^{d + A- a f 



and the mixed part S£u depends both on lattice gauge field and the constrained gauge field. 

J__ r , „ , „ 1 

g 2 a 2 



2* = ^2Tr[d,U r (x)d^(x)} + — 2 gA- a J^ K . 



(2.2) 



(2.3) 



(2.4) 



(2.5) 
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Here the link current 

Jti a NK^) = E ^Trl T a [U r (x)i d+ U} (x) + U}(x - ar)i d + U r (x - ar)} 

Substituting back the expression for A~ a from the constraint equation 

{d + ) 2 A -a = ^ {j +a K _ jr)) 



with 



7+ a (x)=2T7 T (x)r a 77(x 



(2.6) 



(2.7) 



(2.8) 



where 77 is the dimensionless two-component lattice fermion field, in the A a dependent terms in the Lagrangian density, namely, 



2« 2A J * + 8 ( > + 2a 2 Jum 



we generate the terms 



1 



1 



Collecting all the terms, the canonical Lagrangian density for transverse lattice QCD is 



1 



if = ^ 77 t( x )^-T,(x) + J^7'r[^^(x)^C7;(x)]-^-T,t(x)^ T T,(x) 



+im \ T) f (x) O s - \u s (x) tj (x + as) - T7 (x) 



1 r 



77 f (x + ar) U} (x) - T7 f (x) 



- 1 1 r\ 



- -7 f (x + af )t/ r f (x) - tj f (x)] o r ^-a s [U s (x)r] (x + as) - 17 (x)] 
H — ^ V < 7 , r[f/ r (x)f/ s (x + af)f/_ r (x + flf + fls)I/_ 1 -(x + as) - 1 

_ 8 2 j+a ( 1 "\ f +a 

Here ffi = a2, 62 = — G\. In the two-component representation 1 10], the dynamical fermion field 

1 , 



where tj is the dimensionless two component lattice fermion field 



±T7(* ,x x ) 




(2.9) 



(2.10) 



(2.11) 



(2.12) 



2. Effective Hamiltonian 

Because of the nonlinear constraints U^U = 1, det U = 1, it is highly nontrivial to quantize the system. Hence Bardeen and 
Pearson 1 1 ] and Bardeen, Pearson, and Rabinovici 1 2] proposed to replace the nonlinear variables U by linear variables M where 
M belongs to GL(N,'^'), i.e., we replace ~U r (x) — > M r (x). Once we replace U by M, many more terms are allowed in the 
Hamiltonian. Thus one needs to add an effective potential V e ff to the Lagrangian density 



V eff = -^j Tr(M Jf M)+X\ 7>[(M f M) 2 ] + X 2 [det M + H.c] - 



(2.13) 



Thus, the effective Hamiltonian for QCD on the transverse lattice becomes 

P fb = P ffree+ P V + P fhf + P hf + P chnf 



The free fermion part is 



The effective potential part is 



r f free 



+P~ +P~ +P~ +P~ 
qqc i 1 ggc ' 1 qgc ' 1 p • 



J £ /x-^(m 2 + ^)77 t (x)^7 7 ( 



Py = Jdx'a 2 ^ Tr(M' i M)-Xi Tr[(M f M) 2 } - X 2 [det M + H.c}- 



The free helicity-flip part is 



/I 1 

Helicity flip associated with the fermion hop is 

/j j - 
dx~ £ 51 H f ( x ) & - [ Ms W 77 ( x + fl§ ) 



i'mg / dx 52^2 ^(x + ar) Af*(x) 

J vr 



,11 , . 
a d + 



Canonical helicity non-flip terms are 



chnf 



-4 / ^- fl 2 ^^[7 7 t (x + fl f)Mj(x)]a r ^ T a s [77(x)] 

Q J x f v 

-4 / dx-a 2 Y,Y,^( x )^^ T a s [M s {x)ri(x + as)} 
a j x rs , ;c 

s 2 f 1 

/ rf*- fl 2 ££[i7 + (x + a?)Mt(x)]di.^ T ffv[M v (x)77(x + fls)] 



The four-fermion instantaneous term is 



" 2 fj / ^I^ t (x)r77(x)^ I 7? t (x)r''77(x 



The four link instantaneous term is 



The fermion - link instantaneous term is 



P w = 2 8 ^Jdx-^J+ a NK (x)-^r ] \x)T a r ] (x), 



The plaquette term is 



P p = / dx a 2 Y,Y. | Tr M r (x)M s (x + ar)M- r (x + ar + as)M- s (x + as) 



Here 



j link( x ) = H Tr { T a [M r (x)i d + M}(x) +M}(x - ar)i d + M r (x - ar) 
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3. Violations of hypercubic symmetry 

The canonical helicity non-flip interactions given in Eq. ( 12.191 1 for r ^ s break the hypercubic symmetry on the transverse 
lattice. For interacting theory this is also true for the Hamiltonian with symmetric derivative. In the free field limit they do not 
survive for Hamiltonian with symmetric derivative but for forward-backward derivative they survive. In that case, in the free 
field limit they reduce to 



1 



3 /*-LI 



r]\x + ar)d r a s ^-r]{x) 



-77 t (x)d- r di 3 -77(x-|-as) 



- tj f (x + at ) a, a s — tj (x + as) 



Going to the transverse momentum space via 



we get 



T](x-,x l )=Jd 2 k L e ikL - x± A .±CO 



r dx J d 2 k^^{x )a 3 ^ T ^> k ±(x ) 



sin(k y a) — sin(k x a) + sm(k x a — k y a] 



(2.25) 



(2.26) 



(2.27) 



Thus the violations of hypercubic symmetry are of the order of the lattice spacing a. Sign in front of this term changes if we 
switch forward and backward derivatives. 

In our numerical studies presented in Ref. | if] and in this work, we have set the coefficients of hypercubic symmetry violating 
terms to zero. 

B. Canonical transverse lattice QCD with the Wilson term 

When one uses symmetric derivatives for the fermion fields, doublers arise as a result of the decoupling of even and odd lattice 
sites. To remove the doublers one may use the Wilson fermions |7] or the Kogut-Susskind fermions I8l ll ill . In this subsection, 
the details of the structure of the Hamiltonian resulting with the modification of the Wilson term are presented. 



1. Constraint equation 



The symmetric derivative is defined by 



D^ix) = — [f/,.(x)v/ ± (x + af)-f/_ r (x)v/ ± (x-af)]. 
2a 

Again we make the replacement -,U r (x) = M r (x). Anticipating doublers, we can add a "Wilson term": 



(2.28) 



8^ = -\jf(x)[gM r (x)\(f(x + ar)-2\jf{x)+gM^ r (x)\jf{x-ar)} (2.29) 

where K is the dimensionless Wilson parameter. Explicitly, in terms of the dynamical field y/ + and the constrained y/~ 

85? = -V/" t (x)y°[^M,-(x)v/ + (x + af)-2v/ + (x)+gM_ r (x)v/ + (x-af)] 
a 

+ f (x) y° [gM r (x) y- (x + af ) - 2 yr (x) + gM_, (x) yr (x - af )] . (2.30) 



The constraint equation for y/ in the presence of the Wilson term is 

/<9 + y/~(x) = mf\i/ + (x) 
a r 



[gM r (x) y/+ (x + ar) - gM- r (x) y/+ (x - of)] 

2a 



7° [gM r (x) y/+ (x + ar) - 2 y/+ (x) + gM_ r (x) y/+ (x -at)]. 



2. Hamiltonian: Symmetric derivatives and the Wilson term 



After a great deal of algebra, we arrive at the Hamiltonian, 



P f free + P V + P hf + P whf 

+ P clmf + °wn/l + P wnf2 

+P~ +P~ +P~ +P~ 
' 1 qqc ' 1 ggc ' 1 qgc ~ 1 p • 



The free fermion part is 

The helicity flip part is 

P hf = s 



P ffree = Jdx a 2 £l(m + 4^)V(* ,x)^77(^ ,*)• 
/ ^"L| (»J + 4^) ^7? t (x)^d r ^ T [M r (x)77(x + fl f)-M_ r (x)77(x- fl f)] 

7C\ 1 1 

+ 4-J — ^[77 t (x-af)d r M,.(x-flf)-77 t (x + flf)d r M_ r (x + flf)] 



The Wilson term induced helicity flip part 



P whf = g 2 jdx ^|^^^[77 t (x- fl f)M r (x- fl f) + 7 7 t (x + flf)M_ r (x + flf) 



1 

K 1 

a 2a 1 



M s (x) 77 (x + as) - M- s (x) 77 (x - as) 



rj\x — at)a r M r (x — at) - rj Jf (x + ar)a r M_ r (x + ar) 



1 



M s (x) 77 (x + as) + M- s (x) 77 (x - as) 
The canonical helicity non-flip term arising from fermion constraint is 



chnf 



-g 2 I dx YsJ-lHH rf{x-ar)6 r M r {x-ar)-rf{x + ar)6 r M- r {x + ar) 
M s (x) 77 (x + as) - M- s (x)tj (x - as) 



1 

id+ 



The Wilson term induced helicity non flip terms are 



wnfl 



-g j dx £ j (m + 4^j ^77 t (x)^^[M r (x)77(x + fl f)+M_,.(x)77(x-flf)] 
+ (m + 4-J [77 t (x-af)M r (x-af) + 77 t (x + af)M_,.(x + af)] tj+TM > 
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and 



wn/2 



g 2 I dx '? t (x-flf)M r (x-af) + 77 t (x + af)M_ r (x + flf) 



1 



M s (x) 77 (x + as) + M_ s (x) 77 (x - as) 



(2.38) 



Comparing the Hamiltonians with a) forward-backward derivative and b) symmetric derivative with the Wilson term we notice 
that the only differences are in the particle number changing interactions, namely, helicity flip and helicity non-flip terms. 



III. ONE LINK APPROXIMATION 



A. Relevant interactions 



In one link approximation, for either Hamiltonian, the four link instantaneous term and the plaquette term do not contribute 
and only the link mass term of the effective potential contributes. Further, in the case of the forward-backward Hamiltonian, the 
helicity non-flip part proportional to g 2 does not contribute. For the Wilson term modified Hamiltonian, the Wilson term induced 
helicity flip part P~ h p the canonical helicity non-flip term P~ nh f and the term proportional to K 2 in the Wilson term induced 
helicity non-flip part do not contribute. Thus in the case of the Wilson term modified Hamiltonian the entire fermion hopping 
with no helicity flip arises from the Wilson term. 



B. Comparison with one gluon exchange in the continuum 



It is interesting to compare the one link approximation on the transverse lattice with the one gluon exchange approximation 
in the continuum. In the latter, a major source of singularity is the term in the quark - gluon vertex where fe x (k + ) is the 
gluon transverse (longitudinal) momentum. This originates from the A~/+ interaction term in the Hamiltonian via ^r<5 x ■ 
A x contribution to the constrained field A~ . This term gives rise to quadratic ultraviolet divergence in the transverse plane 
accompanied by linear divergence in the longitudinal direction in fermion self energy. On the transverse lattice, d + A~ °= j+J^ INK 
so that A~Jy — ► Jij NK rg+ji Jq- Thus a term which gives rise to severe divergence structure in the continuum gets buried in the 
fermion-link instantaneous interaction term which gives rise to a term in the gauge boson fermion vertex in the continuum in 
Abelian theory. In the non-Abelian gauge theory this gives rise to a term in the quark-gluon vertex and also to the instantaneous 
quark-gluon interaction in the continuum. 

The transfer of the troublesome term from quark-gluon vertex in the continuum theory to quark - link instantaneous interaction 
term in the lattice theory has an interesting consequence. In the continuum theory, the addition of a gluon mass term by hand 
spoils the cancellation of the light front singularity between one gluon exchange and the instantaneous four - fermion interaction. 
On the transverse lattice this cancellation is absent anyway with or without a link mass term. 



C. Longitudinal dynamics and effects of transverse hopping 



We first consider the dynamics in the absence of any link. In this case, fermions cannot hop, and at each transverse location 
we have (1+1) dimensional light front QCD which reduces to the 't Hooft model in the large N c limit. In this case quark and 
antiquark at the same transverse position interact via the spin independent instantaneous interaction which, in the non-relativistic 
limit reduces to the linear potential in the longitudinal direction. The only parameters in the theory are the dimensionless fermion 
mass nif — am and the gauge coupling g. The spectrum consists of a ground state and a tower of excited states corresponding to 
the excitations of the qq pair. 

Next consider the inclusion of the qq link states. There are four independent amplitudes corresponding to whether the quark 
is left, right, above, or below the antiquark. With non-zero mass of the link, these states lie above the ground state of pure quark 
- antiquark system. Further the q, q and link (which are frozen at their transverse positions) undergo fermion - link instantaneous 
interactions in the longitudinal direction which further increases the mass of qq link states. Now the quark or antiquark can hop 
via helicity flip or helicity non-flip. Here we find a major difference between the Hamiltonians resulting from forward-backward 
derivative and symmetric derivative. Let us first consider the helicity flip hopping term in the forward-backward case 



P hf = -imgjdx T7 t (x)d r - — T](x + flf) + T7 t (x + af)M r t (x)d r - — T](x) 



(3.1) 
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If we consider transition from two particle to three particle state by a quark hop, then the first term in Eq. J3 - 1 1 corresponds to 
| 2) — >| 3a) and the second term corresponds to | 2) — >| 3b). The helicity flip term in symmetric derivative case, after making 
some shifts in lattice points, can be written as 



p ;> - -»(- +4 35/* r ?x 



?7 f (x) 6 r (x) 77 (x + at) - 77 f (x) a r M r (x) 77 (x + at) j 

77 t (x)a r ^ r M_,.(x)77(x-flf) - 77 t (x)d- r M_ r (x)^ r 77(x-af) 



(3.2) 



For the Hamiltonian with symmetric derivative, a quark or antiquark hopping accompanied by helicity flip has opposite signs 
for forward and backward hops. On the other hand, hopping accompanied by helicity non-flip have the same signs. As a result, 
there is no interference between helicity flip and helicity non-flip interactions |7]. In the case of the Hamiltonian with forward- 
backward derivative, quark or antiquark hopping accompanied by helicity flip has the same sign for forward and backward hops. 
As a consequence the helicity non-flip hop can interfere with the helicity flip hop. This has immediate consequences for the 
spectrum. In the case with symmetric derivative, in lowest order perturbation theory, the helicity zero states mix with each other 
which causes a splitting in their eigenvalues resulting in the singlet state lower than the triplet state. On the other hand, helicity 
plus or minus one states do not mix with each other or with helicity zero states resulting in a two fold degeneracy. In the case 
with forward and backward derivatives all helicity states mix with each other resulting in the complete absence of degeneracy. 

IV. SINGULARITIES, DIVERGENCE AND COUNTERTERMS 

Since the transverse lattice serves as an ultraviolet regulator, we need to worry about only light front longitudinal momentum 
singularities. 

A. Tree level 



We take all the terms in the Hamiltonian to be normal ordered. At tree level this leaves us with singular factors of the form 
in the normal ordered four fermion and fermion link instantaneous interactions. The singularities are removed by adding the 
counterterms used in the previous work |4] on transverse lattice. The explicit forms of the counterterms are given in Appendix 
IClin the appropriate places. 



B. Self energy corrections 



In the one link approximation, a quark can make a forward (backward) hop followed by a backward (forward) hop resulting 
in self energy corrections. In a single hop, helicity flip or non-flip can occur. In the case of symmetric derivatives, helicity flip 
cannot interfere with helicity non-flip, and as a consequence, self energy corrections are diagonal in helicity space. In the case of 
forward and backward derivatives, the interference is nonzero resulting in self energy corrections, both diagonal and off-diagonal 
in the helicity space. Similar self energy corrections are generated for an antiquark also. These self energy corrections contain a 
logarithmic light front infrared divergence which must be removed by counterterms. In Appendix|E|we present the explicit form 
of counterterms in the two cases separately. In previous works on one link approximation |4, 5, 7], these counterterms were not 
implemented. 



V. NUMERICAL RESULTS 

We diagonalize the dimensionless matrix a 2 P~ . We further divide the matrix elements by g 2 Cf which is the strength of the 
matrix elements for four fermion and fermion - link instantaneous interactions. Now, define the constant G with dimension of 
mass by G 2 = $^C f . DLCQ yields M 2 /G 2 . 

The dimensionless couplings are introduced |4] as follows. Fermion mass nif = m/G, link mass ji/, = ji/G, particle number 
conserving helicity flip coupling m,f/(aG) = ntfCi, particle number non-conserving helicity flip y/Ngm/ / '(aG) — nifC^, and 
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particle number non-conserving helicity non- flip VNg/(a 2 G 2 ) = C 3 . In the case of the Wilson term modified Hamiltonian, 
we have fermion mass term rtif = (m + 4-K/a)/G, helicity-flip coupling \/Ngnif/ (2aG) — tntCi, and helicity non-flip coupling 
s/Ngm/K/ (aG) = irifC^. 

All the results presented here were obtained on a small cluster of computers using the Many Fermion Dynamics (MFD) code 
1 12] implementing the Lanczos diagonalization mathod in parallel environment. For low K values, the results were checked 
against an independent code running on a single processor. 

A. Cancellation of divergences 

As we already mentioned, we encounter jjpry singularities with instantaneous four fermion and instantaneous fermion - link 
interactions which give rise to linear divergences. We remove the divergences by adding appropriately chosen counterterms. 
We have numerically checked the removal of linear divergence by counterterms in DLCQ. First we consider only qq states with 
instantaneous interaction. We study the ground state eigenvalue as a function of K with and without the counterterm. Results are 
presented in Fig. ^( a )- Next we consider only qq link states with fermion-link instantaneous interaction with and without the 
counterterms. The behavior of ground state eigenvalue as a function of K is presented in Fig. [D(b). In both cases, it is evident 
that the counterterms are efficient in removing the divergence. 

B. qq at the same transverse location 

Next we study the spectrum of the Hamiltonian in the absence of any links. Since, in this case, the Hamiltonian depends 
only on the dimensionless ratio -±- we fix g = 1 and vary mj to study the spectra. The Hamiltonian matrix is diagonalized for 
various values of K. The convergence of the ground state eigenvalue as a function of K is presented in Table [I] The ground 
state wavefunction squared as a function of the longitudinal momentum fraction x is plotted in Fig. [2] The convergence of the 
wavefunction has a very different behavior as a function of fermion mass ntf. As can be seen from this figure, the convergence 
in K is from above for heavy w/ and from below for light mj. As a consequence the wavefunction is almost independent of K 
when nif is of order g. 

C. Results of the one link approximation 

We encountered logarithmic infrared divergences due to self energy corrections and, in AppendixlEl we discuss the associated 
counterterms. In Fig. [5]we show the effect of self energy counterterms on the ground state energy in the two Hamiltonian cases 
we studied. 

The quark distribution function for the ground state and the fifth state for the set of parameters m/ = 0.3, fib = 0.2, C2 — 
0.4, C3 = 0.01 and K = 30 is presented in Fig. [4] In this figure we also present separately the contribution from two particle and 
three particle states. As expected, the contribution from the three particle state peaks at smaller x compared to the two particle 
state. The exact location of this peak depends on the link mass. The convergence of lowest four eigenvalues with K for the 
Hamiltonian with forward-backward and symmetric lattice derivatives is shown in Table HP for mf = 0.3, fit, — 0.2. We also 
show the results extrapolated to K — > 00. 

It is interesting to see the effect of fermion - link instantaneous interaction on the low lying eigenvalues. In its absence, there 
is no confining interaction in the longitudinal direction in the qq link sector. Furthermore, the mass of the lowest state in this 
sector corresponds to the threshold mass in this sector. Since its mass is lowered, it mixes more strongly with the qq sector in the 
ground state. The fifth state now corresponds to an almost free qq link state with infinitesimal qq component as shown in Fig. [5] 

VI. SUMMARY, DISCUSSION AND CONCLUSIONS 

In this work we have performed an investigation of qq states using two different light front Hamiltonians in the one link 
approximation. The Hamiltonians correspond to two different ways of formulating fermions on the transverse lattice, namely, 
(a) forward and backward derivatives for and respectively or vice versa and (b) symmetric derivatives for both 
and In the latter, fermion doubling is present which is removed by an addition of the Wilson term. In this case there 
is no interference between helicity flip hop and helicity non-flip hop and, as a result, the qq component of the ground state 
wavefunction which has helicity plus or minus one are degenerate. In the former case, interference between helicity flip and 
helicity non-flip leads to the absence of degeneracy in the low lying spectra. One can recover approximate degeneracy of helicity 
plus or minus one components only by keeping the strength of the helicity non-flip hopping very small. In the case of forward 
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and backward derivatives, terms are also present which violate hypercubic symmetry on the transverse lattice. They become 
irrelevant in the continuum limit when the linear variables M are replaced by non-linear variables U. We have removed them 
entirely from the Hamiltonian in the present investigation. 

Since the one link approximation is very crude, we have not attempted a detailed fit to low lying states in the meson sector. 
Instead, we have explored the effects of various coupling strengths on the low lying spectra and associated wavefunctions. In our 
work, longitudinal dynamics is handled by DLCQ. We have performed a detailed study of various convergence issues in DLCQ 
using a wide range of K values. 

We summarize our results as follows. We have shown the effectiveness of appropriate counterterms in the qq and qq link 
sector to regulate the instantaneous fermion and fermion - link interactions respectively. We have also checked the cancellation 
of logarithmic divergences due to self energy effects. In the limit where fermions are frozen on the transverse lattice but undergo 
instantaneous longitudinal interaction, we have studied the convergence of ground state wavefunction with respect to K for three 
typical values of the fermion mass. We have studied how the presence or absence of fermion - link instantaneous interaction 
in the qq link sector affects the wavefunction of low lying states. We have also studied the consequences of the interference of 
helicity flip and helicity non-flip hopping in the Hamiltonian with forward-backward derivatives. This interference is absent in 
the symmetric derivative case. 

For future studies, we would like to address the problem of mesons containing one light and one heavy quark in the context 
of heavy quark effective theory on the transverse lattice. A systematic study of the effects of sea quarks and additional links on 
the meson observables also need to be undertaken. A major unsettled issue in the transverse lattice formulation is the continuum 
limit of the theory when nonlinear link variables are replaced by link variables. It will be interesting to investigate the light front 
quantization problem with non-linear constraints. In this respect the study of non-linear sigma model on the light front appears 
worthwhile. 
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APPENDIX A: STRUCTURE OF TERMS IN DLCQ 



We use DLCQ for the longitudinal dimension (— L <x < +L) and implement anti periodic boundary condition for the two 
component fermion field, 

rj c (x-,x) = 4fl^ E [fc c (m,x,A)e-^"/( 2L )+4(m,x,-A)e''^"/( 2L )] (Al) 

V2L x m=l,3,5,... 

with 

{fc e (m,x,A),^(m',x',A')} = {d c (m,x,X),dl(m ,x ,X')} = 8 mm/ 8 XjX i8 c , c >8 x ^. (A2) 
The link field has periodic boundary condition (with the omission of the zero momentum mode), 

M rpq (x-,x) = -±= £ -i[B_ rw (m,x+flf)r ta "/ i +fit M (m 1 x) e i -"/ i '] (A3) 

with 

[B r pq (m,x),Bl ts (m',x')] = 8 mm > 8 XtX i 8 r y 8 ps 8 q , . (A4) 

The Hamiltonian P = ^H. 

In the following subsection we give the explicit structure of terms in the Hamiltonian in the forward-backward case in DLCQ 
restricting to those relevant for the one link approximation. 



1. Mass terms 



Mass terms: 



H f free = ™ 2 £EEE - bl(n,x,X)b c (n,x, A) + 4("> x > X)d c (n,x,X) 



x c A, n 



x r n 



B T r (m, x)B r (m,x) + B^_ r (m, x + ar)B- r (m, x + af) 



2. Four fermion instantaneous term 



The four fermion instantaneous term which gives rise to a linear potential in the color singlet state 

g 2 

E E E^a'Sw E 

;t " cl j c"c"'XX'X"X'" x mim2»'3'n4 

x bl(mi,x,k)dl m (m4,x,—X"')b l j(m2,x,k')d c ii(m3,x,X'") 



(m 3 -m 4 ) 2 mi+m4 ' m2+m3 • 



3. Helicity flip terms 



Particle number conserving terms: 



" r x A^Aa mi " n 

[bl(m l ,x,Xi)b c (m u x,X 2 )+de(m u x,-X2)d c (m u x 7 -X 2 )] . 
Particle number non conserving terms: a typical term is 

mg 1 VV V t - V — 1 s 



a r x XiM 



m\tn2tnT, 



m3 2m3 + m2 



frj (mi , x, X\ )B_ rcc i (mo, , x + ar)b c i (ni2,x + ar, X 2 )- 



4. Helicity non flip terms 



Two operators: 



2 1 

~tEEE~ [bi{n,x,X)b c (n,x,X)+d c '(n,x,X)dc{n,x,X] 



"TIT" 



Three operators: 
A typical term is 



11 11 

:EEE E - 



b c (mi,x,X)B_ rcc i (m3,x + ar)b c i (w2,x + af,A). 
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5. Fermion - link instantaneous term 

A typical term is 

2 4^3EE E YJcd T M' E ~i^.~Jwa 

* %a x r cdd'dd! mi m 2 m 3 m 4 V m 3 V m 4 

b\{m x ,x, X\ )b d i (m 2 , x, A 2 )B_ rc ' c " (w 3 , x + af )si rc „ c (m 4 , x + af ) 

(-)(m 3 +m 4 )/(mi-m 2 ) (All) 

APPENDIX B: STATES IN DLCQ 

We will consider states of zero transverse momentum. In the one - link approximation, the gauge invariant states are qq states 

V^V VV d y(9) m 

bfa u j(q),Oi) dl{n 2 ,y{q),a 2 ) | 0) (Bl) 



and the link states 



= N ~7T7 ~7! E E E E E 5 y('),yW 5 y(<?),y(?)-«s 
« V2 dd , s y{g) m y(;) 

bl{m,y{q),Oi) B] dd ,{m,y{l)) dl,(n 2 ,y(q),o 2 ) | 0) 



and 



1 1 1 



= W ^5 E E E E E 5 y('),y(9) 5 y(q),y(q)+as 
v VK vz rfrf' s y(g) y(g) y(0 

ftj(ni,y( 9 ),ai) B^ rf ,(« 3 ,y(0) 4("2,y(?),^) I 0). (B2) 
We shall consider transition from these initial states to the following final states: The qq state 

(0 I d e {n' 2 ,z{q),a' 2 ) b e {n[,z{q),a{) (B3) 



and the qq link states 



^ ' = ^ ~7i7 ~7?E E E E E 5 </),<9) 5 z( ? ),z(?)-ai 

(0 I d e (i^,z(§),<^) B tee ,{n' 3 ,z(l)) M»i>*(«W) ( B4 ) 



and 



1 1 1 



Vf Ml ee > t z ( q ) z (q) z (/) 

(0 | d e (ri 2 , z{q) , a 2 ) B_ tee , (n' 3 , z{l) ) b^ (n[ , z{q) , a[) 



(B5) 
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APPENDIX C: FORWARD-BACKWARD DERIVATIVES: MATRIX ELEMENTS IN DLCQ 
1. Transitions from two particle state 

a. To two particle state 
Let us consider transitions to the two particle state: We have, from the free particle term, 

(2' | H f free \2)=m 2 (— + -) Jf 2 (CI) 

where 

^ = <?»!,»'! K,n' 2 8 C 2 ,^ ■ ^ 

From the four fermion instantaneous term we get 

„2 ! 



(2' | H qqc | 2) - 2^ 2 C f 8 ni+mA+n , 2 ^ _ b , )2 



where C 



iv 2 -i 



/ — 2N 



To implement the regulator prescription for jp^i, we add the counterterm matrix elements 

„2 j 



(2' I #cr | 2) = 2^ C/ S niWi+n , J£ ^ (ni _„ /oop)2 5oi,oi W (C4) 



Here the term «/ oo; , = «i is dropped from the sum. 
From the helicity flip term we get 



(2' | H hfl I 2) = - 2 ^L ^ 4; & Zoi ^ + ^ z! a2 a s ^ 5^; 



with 



From the helicity non-flip term we get 



(C5) 

•A/ = 5 «l,n' 1 S n 2 ,n' 2 ■ ( C6 ) 



(2' | //, n/ (l) | 2) = 2^ f- + - ) ^2 • (C7) 
a 1 \n\ n 2 



To three particle state 
c. To the state \ 3a) 



From the helicity flip term we get 



°n2,n> 2 -, 7=r 2- 2- d z(</),y(<?)-at 

"l V»3 «(,) y( 9 ) 

. mg r- 1 1 1 f „ 



14 



From the helicity non-flip term we get 

(3a' | H hnf {2) | 2) 



1 m 1 1 1 st st 



<*n[+2ri 3 ,ni 1 



f I E I 5 z( 9 ),y( 9 )-at 



"l V n 3 » «(,) y( 9 ) 



/jE E E 5 z(?),y( 9 -)+at 



(C9) 



From the helicity flip term we get 

(3&' | // V2 I 2) 



From helicity non-flip term we get 

(3b' | H hnf (3) | 2) 



d. To the state \ 3b) 



m S fT, 1 1 1 V t - x 



"n 2 ,ri 2 



"3 z( 9 ) y( 9 ) 



z(«)j(?)+«t 



a V V2 V47T ^ 



^+2h 3 ,h 2 1 



"1 '"I 



f I 1 5 > 



"3 y(q) 



z(§),y(?)-at ■ 



"1 



1 / — 11 1 



"3 r z( 9 ) y( 9 ) 

8„ 



z(</),y(<?)+«t 



v V2 a2 ^ 2 <Tl ' CTl 

fl E E ^ „ 



(?),y(?)-at 



(CIO) 



(Cll) 



2. Transitions from three particle (a a link) state | 3a) 

a. To three particle state 



From the free particle term, we get 



(3a'\H free \3a)= ( m 2 (- + -) + l - l i 2 - W 3 



with 



•^3 - 5 nij „/ 5„ 2 ^ <5„ 3 y 3 5 CTl <T / S^ o-. 



'2 



(C12) 



(CI 3) 



Diagonal contribution from the four fermion instantaneous term to the three particle state vanishes due to the vanishing trace 
of the generators of SU (N). 
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Contribution from the fermion - link instantaneous term 



8 2 1 

(3a | H qgc (l) | 3a) = -— -~C f S„ 1+2 / +2 „/ S n 2 ,n' 

71 d 



1 (ni-n / 1 +4w 3 ) J_ 

VHi v /»i-»'i+2«3 («i -»'i) 2 V2 ^ ^ 



.El 

% a 



2^f 8n 2 +2n 3 ,n' 2 +2n' 3 ^m,n\ 



(n 2 -n 2 + 4n 3 ) 1 

; rr^ 7= 0„. „i 0„„ „l . 



Counterterm matrix elements in DLCQ to implement the regulated prescription for 



(C14) 



(3a' | H CT (l) | 3a) = — -^C f <5„ 1+2ll3 y i+2 ^ S„ 2j „/ 8 aua , 8^ 

1 (»1 -»/oop + 4w3) 1 



w 1 ma.Y 



I 

->o p- 

+ I 



„ ;oop= i V^v 7 "! -»/oop + 2 »3 ("l-n/« v ) 2 V2 

nimax 1 (»2-»/ op + 4» 3 ) 1 

1 V«3 ^"2 - W/oop + 2«3 ("2 - n/oop) 2 V2 



n loop- 



where n lma3C < «[ + 2n 3 and n 2mfll < n 2 + 2n 3 . 

The contribution from the helicity flip term that conserves particle number is 



(3a'\H hf (l)\3a) 



? _ £ £ £ 



The contribution from the helicity non-flip term that conserves particle number is 

(3a'\H hnf (l)\3a) = ^(- + -\yV 3 . 

a L \n\ n 2 I 



(C15) 



(C16) 



(C17) 



b. To two particle state 



From the helicity flip term we get 



<V,,Hi+2n 3 1 V V £ 

°n 2 W 2 — L L <fy<?),y(<?)+«s 



From the helicity non-flip term we get 



y(«) 

1 1 1 



_ __ Y t x 



°nun[ ~, 7= L L °z(q),y(q)-aS ' 



"2 V"3 z((? -) y(9 -) 



(2' I tf ftn/ I 3a) = -gl i 4= -i= 4 .,■ 5. 



'a 2 V V2 V4tt 



CTl,CT "(72,^9 



(C18) 
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^Vj,ni+2« 3 1 



m 

1 / — 11 

VN 



EE I* 



W<?)j(<?)+"S 



' 5„ 



0„. „/ j — ^ 2- d z(q),y(q)-as 



" 3 y(?) 



(C19) 



3. Transitions from three particle (q q link) state | 3b) 
a. To three particle state 



From the free particle term, we get 







1 2 M 




•"4 






«2' 


2 «3 / 



with 



5 ni,n'j 5 » 2 ,n' 2 5 «3,* 3 5 CTi,ct( 5 CT 2 ,^ 



(C20) 



(C21) 



The diagonal contribution from the four fermion instantaneous term to the three particle state vanishes due to the vanishing 
trace of the generators of SU(N). 

The contribution from the fermion - link instantaneous term is 



(3b'\H qgc (l)\3b) = -s--c f 8 n 



5? 

11+2/13, II j+2rt4 ^2; n 2 



1 («i-» / 1 +4» 3 ) J_ 8 8 

^ v /»i-»'i+2»3 («i-"'i) 2 V2 ^ ° 2 ^ 



r i 

1 



(«2 — n'2 + 4«3) 1 



5„. „/ 5„ 



>JV"2-"2 + 2"3 (»2-« 2 ) 2 ^ ^ ^ 



(C22) 



Here also we have the counterterm matrix elements given in Eq. (IC15t . 

The contribution from the helicity flip term that conserves particle number is 



(3b'\H hf (l)\3b) = -2-d nun , 8 mA 8 niA 

7" E xl>°rXo { 8 aia > + 7 E X-a 2 VrX-a< ^o' 



«i 7 1 "2 T 

The contribution from the helicity non-flip term that conserves particle number is 



(3b'\H hnf (l)\3b) 



1 1 

— + — 



(C23) 



(C24) 



b. To the two particle state 

From the helicity flip term we get 



(2' \H hf2 I 3.) = f VN 11,^x1, a sXoi 8^, 
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1 — r I & ' 5 / . _ / 

«i n,l "2' n 2 «i+2« 3 ,«i 

/ JC\ JC / — 1 1 1 
- m + 4- )-Vn ^^=5 ,8 i 

E E E 5 z( 9 -)j(9)+«t 

' y(«) 

' 4 + — ) 5 , 5/ , , . (D2) 



fe To ?fre stare | 3b) 



Helicity flip: 



Helicity non-flip: 



(3*115,12) = (,„ + 43l^-i=E^*fe^ 



1 1 1 



E E 5 z( 9 ),y(g)+ai 

y(«) i(«) 

x «1 rtj / " 2 ^2 «l+2« 3 ,«l 
KT\ 1 /-l 1 1 



E E 5 z(?),y(,?)-ai 

y(« 

« 2 « 2 / "''"1 "2+ 2 «3>"2 



(3«'|P- n/ i|2) 



, K\ K /— 1 1 1 

m+4— - vW = ^=<5 / <5 / 

a/ a VV2V47r <Tli<T i 



E E E 5 z( 9 ),y( 9 )+at 

' y(?) z (?) 

1 / 1 1 \ - 



« 3 \ 1 / 

K"\ K" ,— \ 1 1 

m + 4- — VN =^=8 i 8 i 

a) a V V2 V47T ^' ff 2 

E E E 5 z(?),y(5)-at 
' y(«) *(?) 

— ' 4- + — ) 5 ' 5 / ^ / . 

v « 2 «2 / "''"l "2+ 2,, 3."2 



(D3) 



(D4) 



2. Transitions from three particle state | 3a) to two particle state 



Helicity flip: 



« r *w - (-HI) 5^735 



l l l 



E E ^(?),y(9)+a§) 



E E 5 z(«),y(5)-fl§) 
z(<?)y(<?) 



Helicity non-flip: 



(2' | p- nfl | 3a) = - (m + 4-) -VjV-^= -t=5„ » 5 , 
x 1 ' V a/ a V V2 V4tt ffl,<T i 

EEE 

s *(?)y(«) 

1 Z' 1 i A o 

h -r 5 i 8 „ / 

V a) a V V2 V47T 2 i' 1 

E E ^(?),y(«)-hrf) 

y(q) 

— (UA 3,8,,. 



3. Transitions from three particle state | 3b) to two particle state 



Helicity flip: 



Helicity non-flip: 



1 1 1 



a 2 ,a 2 



E E ^(?),y(«)-as) 
z(<?) y(«) 



^"3 Ul «! / " 2 '"2 "1+ 2 «3,« 1 



E E 5 z(q),y(q)+as) 

<q)y(q) 



V"3 \«2 "2/ "2 +2 "3'"2 



1 1 



1 



EEE 5 zW,y(<?)-«§) 
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-== [ — + 4- I <? i 8 , , / 

V"3\"i nJ " 2 '"2 "i+ 2 "3"i 

— V^V— — f= —7=8 i 8 i 
a) a V y/2 CTl>CT i 

III 4(?) >y ( ? )+ a s) 
s *(?) y(?) 



— I - + - 

V"3 \ n 2 " 2 / " ni '"l ~ n 2+ 2n 3-."2 



— (t + -| 5 '5 , . (D8) 



APPENDIX E: SELF ENERGY COUNTERTERMS 

In this Appendix we list the self energy counterterms. 

1. Symmetric derivatives case 

The counterterm for self energy for a quark or an antiquark with longitudinal momentum n\ due to double helicity flip hops 

2 1 («i-«'i) 2 

CT\ — — y — = ~-/ 7TT- (El) 

n\ rty /i z ni«j -\-mr(n\ — n\) z 
The counterterm for self energy for a quark or an antiquark with longitudinal momentum n\ due to double helicity non-flip hops 

CT2 = l£l , l-; + f , v , (E2) 

2. Forward and backward derivative case 

In this case we have three types of contributions: (1) helicity flip acting twice, (2) helicity non-flip acting twice and (3) 
interference of helicity flip and helicity non-flip hops. The first two are diagonal in helicity space but the last one is off-diagonal 
in helicity space. 

The transition from state | 2) to state | 3a) and back due to a quark hop gives rise to longitudinal infrared divergence. In this 
case the counterterm due to double helicity flip is 

1 

Cr 3 = 2 £ -— j— -L (E3) 

The counterterm due to double helicity non-flip is the same without the factor of 2. The transition from state | 2) to state | 3b) 
and back due to a quark hop does not give rise to longitudinal infrared divergence. Similarly the transition from state | 2) to state 
| 3a) and back due to an antiquark hop does not give rise to longitudinal infrared divergence. The transition from state | 2) to state 
| 3b) and back due to an antiquark hop gives rise to longitudinal infrared divergence which requires counterterms the explicit 
forms of which are the same as in the quark case for the transition from | 2) to state | 3a). Lastly we consider counterterms 
for self energy contributions arising from the interference of helicity flip and helicity non-flip hopping. The counterterms have 
the same structure as in the case of helicity non-flip transitions accompanied by the following extra factors. Since we have two 
possibilities namely helicity flip followed by helicity non-flip and vice versa and these two contributions are the same, we get a 
factor of two. We also get a factor G^Xs where s(s') is the initial (final) helicity and d 1 = a 2 , 6 2 = — c 1 . 
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K 


Eigenvalue (^# 2 ) 


nif — 0.3 


m f = 0.9 


nif — 3.0 


10 


0.620 


4.547 


39.233 


18 


0.693 


4.664 


39.861 


30 


0.745 


4.724 


40.053 


50 


0.788 


4.762 


40.163 


78 


0.819 


4.783 


40.220 


98 


0.832 


4.791 


40.241 




0.869 


4.820 


40.285 



TABLE I: Ground state eigenvalue (in units of G 2 ) for qq sitting at the same transverse location. 





Forward-backward 




Symmetric 






(C 2 


= 0.01, C 3 =0.4) 


(C 2 = 0.1, C 3 = 0.4) 


K 




J% 




Ji} 


Jtl 








10 


0.38041 


0.4800 


0.4899 


0.5996 


0.3486 


0.4507 


0.4507 


0.5980 


18 


0.3722 


0.4968 


0.5110 


0.6447 


0.3402 


0.4673 


0.4673 


0.6409 


30 


0.3606 


0.5027 


0.5210 


0.6680 


0.3288 


0.4702 


0.4702 


0.6620 


42 


0.3511 


0.5029 


0.5240 


0.6765 


0.3189 


0.4677 


0.4677 


0.6682 


50 


0.3457 


0.5019 


0.5246 


0.6790 


0.3130 


0.4651 


0.4651 


0.6693 




0.3243 


0.5022 


0.5313 


0.6979 


0.2913 


0.4589 


0.4589 


0.6837 



TABLE II: Lowest four eigenvalues (in units of G 2 ) in one link approximation. 
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FIG. 1: Effect of counterterm on the ground state eigenvalue, (a) With and without the counterterm in the qq sector for my = 0.3. (b) With 
and without the counterterm in the qq link sector for mj = 0.3 and \L h = 0.2. 
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FIG. 2: Quark distribution function | | 2 of the ground state in the qq approximation for three choices of quark masses with coupling 
constant g = 1.0. 
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FIG. 3: Effect of self energy counterterms on the ground state eigenvalue in the case of (a) symmetric derivative with £2 = 0.4, C3 = 0.1 and 
(b) forward-backward derivative with C2 = 0.4, C3 = 0.01. mj = 0.3, = 0.2 for both cases. 




FIG. 4: (a) Quark distribution function | \j/(x) | 2 of the ground state in the one link approximation, (b) qq contribution to the ground state, (c) 
qq link contribution to the ground state, (d) Quark distribution function | y(x) | 2 of the fifth eigenstate in the one link approximation, (e) qq 
contribution to the fifth eigenstate, (f) qq link contribution to the fifth eigenstate. The parameters are my = 0.3, fib = 0.2, C2 = 0.4, C3 = 0.01 
and A: = 30. 




FIG. 5: Without the fermion - link instantaneous interaction: (a) Quark distribution function | \jf(x) \ 2 of the ground state in the one link 
approximation, (b) qq contribution to the ground state, (c) qq link contribution to the ground state, (d) Quark distribution function | \jf(x) | 2 of 
the fifth eigenstate in the one link approximation, (e) qq contribution to the fifth eigenstate multiplied by 10 4 , (f) qq link contribution to the 
fifth eigenstate. Parameters are the same as in Fig. [4] 



